Biomedical Optics Express

High-resolution microscopy techniques are used across research and industry to analyse biological systems, from biomolecules to subcellular organelles, multicellular models and tissues. As multimodal imaging workflows and quantitative analysis of bioimaging data become increasingly widespread, there is a growing need for materials and methods to calibrate imaging systems and evaluate the fidelity of generated image data. Here, we present three-dimensional microscopy phantoms fabricated using two-photon photolithography from transparent resins that exhibit both broadband visible autofluorescence and Raman scattering across the fingerprint and C-H stretching regions. Suitable for analysis using optical profilometry, the phantoms were dimensionally calibrated with SI traceability using a metrological confocal microscope. Immersible in air and common aqueous imaging media, the phantoms are compatible with a wide variety of optical microscopy techniques, including one and two-photon excited fluorescence and coherent Raman scattering microscopy. We employed a forked wedge design to validate image deconvolution results and a stacked lattice phantom to recover image distortion matrices under realistic biological imaging conditions. We demonstrate the impact of correcting chromatic offsets and axial scaling errors for a representative application: analysis of a cell seeded scaffold using confocal laser scanning fluorescence microscopy. These phantoms provide a versatile platform for calibration, quality control and validation of multimodal imaging pipelines and improved quantitative optical microscopy.

Background: Quantitative lipid assessment is central to identifying rupture-prone coronary plaques and represents a therapeutic target for lipid-lowering therapy. Near-infrared spectroscopy (NIRS)-derived lipid core burden index (LCBI) is well validated and widely used for detecting lipid-rich lesions. Optical frequency domain imaging (OFDI) is increasingly adopted for guiding percutaneous coronary intervention (PCI) due to its high-resolution structural imaging capabilities. Depolarization-sensitive OFDI (depOFDI) provides intrinsic lipid contrast and may enable combined structural and compositional plaque characterization within a single OFDI-based platform. Objective: To define an OFDI-derived lipid metric and evaluate its agreement with NIRS-derived LCBI. Methods: Thirty-three patients underwent both polarization-sensitive OFDI and NIRS-intravascular ultrasound imaging during PCI. After exclusion of 4 datasets, 29 co-registered pullbacks were analyzed. A signal-to-noise-corrected depolarization metric was used to identify lipid-rich regions and generate depOFDI chemograms. maxLCBI4mm value and location, as well as total LCBI, were computed and compared with NIRS. Results: depOFDI demonstrated strong agreement with NIRS, showing high correlation for maxLCBI4mm (r^2 = 0.862) and total LCBI (r^2 = 0.867), along with strong spatial concordance for the location of the maxLCBI4mm (r^2 = 0.900). Bland-Altman analysis of LCBI4mm showed minimal bias (10.7) with 95% limits of agreement of [81.4 to 102.8]. Conclusions: depOFDI enables accurate quantification of lipid burden alongside the high-resolution structural information inherently provided by OFDI. Because depolarization metrics can be derived from polarization-diverse detection available in many commercial OFDI systems, this approach provides a practical pathway toward comprehensive plaque characterization within existing PCI workflows, without the need for additional imaging modalities.

1.We report a repetition-controllable gain-managed nonlinear fiber amplifier (GMNA) that delivers near-infrared 50-fs pulses with pulse energies up to 150 nJ and a widely tunable repetition rate from 1-20 MHz, while maintaining stable pulse quality across the full range. Using this source, we demonstrate label-free multiphoton imaging--including metabolic autofluorescence (2PF/3PF), second/third-harmonic generation, and Simultaneous Label-free Autofluorescence Multiharmonic (SLAM) microscopy imaging--across live cells, human lung spheroids, and hard tissues. We further assess the impact of laser repetition rate on photodamage at fixed pulse energy, supported by preliminary measurements indicating lower damage at lower repetition rate. Collectively, the compact architecture and repetition-rate agility of the GMNA enable real-time optimization of imaging speed, depth, and sample safety for advanced biological microscopy.

Light sheet fluorescence microscopy enables volumetric imaging with high imaging speed, optical sectioning capability, and reduced photobleaching and phototoxicity, and has become a workhorse in bioimaging. However, widely adopted Gaussian light sheets face an inherent trade-off between axial resolution and field-of-view due to diffraction. State-of-the-art nondiffracting light sheets--including Bessel beam, Airy beam, and lattice light sheet--alleviate this trade-off but suffer from optical aberrations that compromise performance with increasing imaging depth. While the integration of adaptive optics offers a promising solution, such integrated systems are typically complex, expensive, and slow due to the need for serial mapping and correction of spatially varying aberrations across the specimen. Here, we present polarization-engineered aberration-resilient light sheet (PEARLS), a new class of monochromatic nondiffracting light sheet with temporally invariant profile and robustness to optical aberrations. In comparison with existing light sheets, PEARLS showed significantly reduced photobleaching and enhanced aberration-resilience, permitting imaging of three-dimensional subcellular dynamics in optically complex environments. We applied PEARLS for noninvasive observations of biological dynamics in various living systems, revealing phenotypic diversity across spatial and temporal scales--from rapid membrane dynamics and organelle interactions in cultured cells to coordinated mitosis and cell migrations in developing embryos.

The spatial resolution of optical imaging systems is fundamentally restricted by the diffraction limit. However, in widefield live-cell microscopy, the achievable resolution is further constrained by the specimen motion, which indicates the existence of a fundamental spatio-temporal resolution trade-off between signal accumulation during the full frame integration and the resulting motion blur. To improve the fidelity with which moving objects can be imaged, a quantitative understanding of this spatio-temporal trade-off is necessary. Here, we present a systematic analysis of motion-induced resolution dynamics measured with spectral signal-to-noise ratio (SSNR). We developed a simulation framework which models the image formation of objects undergoing arbitrary motion, to evaluate the degradation of the spatial resolution under translational and rotational dynamics. Our results demonstrate that for translating objects, the spatial resolution is anisotropically reduced as a function of the orientation of the object relative to the motion vector, leading to the spectral signal-to-noise ratio degrading by up to 50% and the resolution by up to 40% for a 90{degrees} change in the motion direction. Furthermore, we show that for rotational motion, conventional radially averaged metrics such as the Fourier Ring Correlation are not able to quantify the effects of angular blur. On the other hand, the SSNR is able to accurately quantify this degradation. These findings underscore the necessity of an object-oriented imaging approach, in which acquisition parameters such as exposure time are tuned to specific biological spatio-temporal characteristics to optimize the trade-off between motion blur and spatial fidelity.

Alterations of fluorescence properties in retinal pigment epithelium (RPE) cells caused by diseases such as age-related macular degeneration (AMD) highlight the need for detailed analysis of the fluorescent RPE granules at the individual level. Precise segmentation and classification of these granules remain challenging due to their limited visual separability. In this study, we present Classi4RPE, a computational algorithm designed to accurately segment RPE granules and classify them into three categories -- lipofuscin (L), melanolipofuscin (ML), and melanin (M) -- based on fluorescence lifetime imaging data, which provide distinctive contrast. The method is implemented in a custom Python framework and employs seeded watershed segmentation to isolate individual granules. Lipofuscin granules are identified as hyperfluorescent structures with longer lifetimes, while granules with shorter lifetimes are further analyzed based on their spatial lifetime distribution from the center to edge, enabling discrimination of ML from other melanin-rich granules. Our approach achieves high performance, with mean sensitivities of 0.99 for L granules and 0.90 for ML granules, and corresponding specificities of 0.93 and 0.98, respectively, compared to manually annotated ground truth. These results demonstrate the potential of Classi4RPE to surpass human visual limitations and provide a robust tool for quantitative RPE analysis.

Coherent structured illumination microscopy (c-SIM) is a synthetic aperture optical technique for sub-diffraction limit imaging that extends the utility of traditional SIM to non-fluorescent samples. Here, we present a complementary 5-beam implementation of c-SIM that provides enhanced optical sectioning compared to conventional quadrupolar illumination. Since our approach detects intensity images due to coherent light scattering, it avoids the complications associated with detecting complex fields. Through comparative measurements on calibration samples and live microalgae, we show that 5-beam c-SIM effectively suppresses coherent defocus effects, improving image quality while simultaneously providing a 2-fold lateral resolution improvement.

Stimulated emission depletion (STED) microscopy is a super-resolution fluorescence imaging technique that achieves high spatial and temporal resolution by exploiting stimulated emission to induce fluorescence depletion (FD) and is expected to have substantial utility for imaging applications using fluorescent proteins. However, the compatibility of fluorescent proteins with STED microscopy systems has been understood primarily through empirical observations, and there is no established methodology for the rational selection of fluorescent proteins for STED microscopy. In this study, we systematically evaluated the compatibility of commonly used fluorescent proteins with STED microscopy systems by measuring FD properties using transient absorption spectroscopy and fluorescence dip spectroscopy, both of which are classified as two-color spectroscopy (TCS). Fluorescent proteins identified as compatible with the STED microscopy system based on the TCS measurements were employed for three-dimensional STED imaging of cellular samples expressing each protein. In all samples, three-dimensional spatial resolution was improved relative to confocal laser microscopy, with particularly marked improvements in z-axis resolution. These findings demonstrate that measurements of FD properties via TCS provide a robust approach for evaluating the compatibility of fluorescent proteins with the STED microscopy system and for selecting suitable fluorescent proteins for STED imaging.

Inverse-scattering methods enable label-free, quantitative visualization of a samples three-dimensional (3D) refractive index (RI), providing intrinsic and volumetric morphological contrast without exogenous labels. This is achieved by developing computational frameworks that reconstruct the samples 3D RI from a series of scattering measurements acquired under different data-capture conditions. Recent advances have demonstrated successful 3D RI reconstructions in multiple-scattering samples using angle-varying illuminations; however, these studies have primarily focused on non-absorptive samples. Here, we extend the multi-slice beam propagation (MSBP) inverse-scattering framework to reconstruct complex-valued RI, encompassing both the samples conventional RI (real part) and absorptivity (imaginary part). We show that reconstructing complex-valued RI makes the inverse problem ill-posed under angle-varying illumination alone, and that incorporating measurement diversity from both angle-varying illumination and sample defocus is necessary to ensure stable and accurate convergence. Experimental demonstrations were conducted on 1) dyed microsphere samples to characterize accuracy of reconstructed RI and absorptivity; and 2) diverse absorptive scattering samples to demonstrate biological utility. These results represent an important step for label-free volumetric imaging in biological tissue, which typically exhibits both scattering and absorption.

Light-sheet fluorescence microscopy (LSFM) has revolutionized biological imaging by enabling high spatial and temporal resolution with minimal photodamage. However, conventional LSFM techniques often suffer from striping artifacts in the resulting images due to light scattering and absorption within samples, leading to uneven illumination that negatively impacts the accuracy of subsequent image analyses. To address this limitation, we introduce dodecagon light-sheet fluorescence microscopy (dodecaLSFM), a novel approach that maximizes angular diversity to achieve homogeneous illumination and suppress striping artifacts. dodecaLSFM employs diffraction optics and cylindrical lenses to generate twelve light sheets, providing 360 degree omnidirectional illumination that significantly enhances illumination uniformity compared to traditional mSPIM, mDSLM, and ultramicroscopy systems, which use only one or two illumination planes. We demonstrate the effectiveness of dodecaLSFM by achieving high-resolution imaging of whole mouse brain vasculature following tissue clearing, allowing precise morphometric analysis of vascular networks without striping artifacts. Furthermore, we show that combining dodecaLSFM with expansion microscopy (ExM) enables whole-organ 3D imaging at cellular resolution. This novel approach provides an advanced, scalable solution for large-volume imaging, facilitating detailed structural and functional studies across diverse biological applications.

Fluorescent voltage-sensitive dyes (VSDs) enable non-invasive, high-throughput optical measurement of membrane potential in living cells, but the analytical reliability of such measurements depends critically on whether the dye and associated imaging conditions perturb the system under study. Here, we systematically characterise the photophysical performance and cell-perturbing effects of FluoVolt, a widely adopted VSD, across cancer cell lines (GIN31 glioblastoma and SK-MEL-30 melanoma) and primary human macrophages. Photobleaching kinetics were strongly cell-type-dependent, with SK-MEL-30 cells exhibiting complete fluorescence loss within 400 seconds under standard widefield conditions. FluoVolt staining combined with laser excitation caused an approximately 2.5-fold increase in cell detachment relative to unstained controls, and dual-wavelength excitation (488 + 405 nm) reduced GIN31 cell viability by approximately 17.5%. Critically, morphological changes, a transition from elongated to amoeboid-like phenotypes, were detected under staining conditions alone, prior to any laser exposure, indicating baseline dye-induced perturbation independent of phototoxicity. Halving dye concentration and loading time significantly attenuated these effects while preserving measurable fluorescence signal. These findings identify FluoVolt staining and excitation as previously uncharacterised sources of systematic measurement artefact and provide practical, actionable guidance for protocol design, control selection, and data interpretation in optical membrane potential imaging.

Minimally invasive surgery is a powerful technique that enables operations deep within the body while minimizing patient trauma and recovery time. Optical endoscopes are key to providing intraoperative vision but still face challenges due to the loss of essential senses, including depth perception and tactile feedback for tissue evaluation. Thus, it is critical to develop endoscopic imaging technologies that can augment operators with critical information. In this work, we explore a prototype multimodal 3D imaging endoscope that integrates volumetric light-field imaging with laser-speckle contrast imaging to simultaneously capture 3D structure and blood-flow information in a clinically relevant form factor.

Dark adaptation is the essential process that restores visual sensitivity following exposure to bright light, yet the underlying mechanisms remain incompletely understood. Here, we propose a method for assessing dark adaptation in cones using optoretinography (ORG) based on adaptive optics optical coherence tomography (AO-OCT). ORG quantifies cone functional response by monitoring nano-scale changes in the cones outer segment occurring over hundreds of milliseconds after visible stimulation. This method consists of sequential measurements of stimulus-evoked cone responses over the course of minutes of dark adaptation. Each response captures optical path length changes in single photoreceptor outer segments over milliseconds during a multi-minute recovery period following a strong photopigment bleach. We parameterized cone ORG responses and proposed an exponential model linking ORG dynamics to pigment regeneration. Parameters of the ORG response exhibited exponential decay behavior during dark adaptation, and were thus fit with exponential functions and quantified by the resulting decay parameter{tau} . Parameters capturing the amplitude of the ORG responses recovered more slowly than those capturing temporal dynamics of the responses. This difference is consistent with distinct contributions from photopigment regeneration and downstream phototransduction processes. Recovery speed varied by two-to threefold among three normal subjects, suggesting substantial inter-subject physiological diversity. Processes within the cone, including pigment regeneration, are thought to underlie the gains in photopic visual sensitivity that occur in the dark. These findings highlight ORG as an objective and sensitive assay of those cellular mechanisms. While the ORG itself has shown promise as a biomarker of the health of the photoreceptor response to light, the results of this study show that it may also be useful for probing the health of the intra- and intercellular homeostatic mechanisms that support it.

Significance: Continuous, noninvasive, and depth-resolved monitoring of blood-flow-related tissue dynamics remains an important unmet need. Speckle-contrast optical spectroscopy (SCOS), including interferometric implementations such as iSCOS, provides a scalable optical route to blood-flow sensing, but conventional continuous-wave approaches lack intrinsic depth selectivity. Time-of-flight (TOF) gating offers a way to separate superficial and deeper dynamic contributions in layered tissues, such as skin-muscle or scalp-cortex, by resolving photon path lengths. Aim: We introduce a swept-source, single-channel implementation of interferometric speckle-contrast optical spectroscopy (iSCOS) to obtain TOF-resolved temporal speckle contrast, {kappa}^2, from the measured field autocorrelation g_1, and evaluate its feasibility for depth-resolved blood-flow sensing. Approach: A swept-source iNIRS system operating at 780 nm acquired interferometric signals, which were Fourier-transformed along the optical-frequency axis to recover complex TOF-resolved speckle fields. Temporal speckle contrast was then estimated at each TOF gate indirectly from g_1 using the speckle-visibility relation. Diffusion-based numerical simulations were first used to compare the direct variance-based estimator and the indirect g_1-based estimator under varying reduced scattering coefficient, diffusion coefficient, additive noise level, and bi-layer geometry. Because the simulations showed that the g_1-derived {kappa}^2 estimator was substantially less sensitive to additive noise than the direct estimator, this estimator was used for the main phantom and in vivo analyses, while the direct estimator served as a simulation comparator. The g_1-derived estimator was then applied to liquid and bi-layer phantoms, followed by proof-of-concept in vivo measurements on the human forearm during cuff occlusion and on the forehead during a Sudoku task. Results: TOF-resolved kappa2 curves recovered with the g_1-derived estimator matched DWS theory across scattering coefficients, photon path lengths, and exposure times. The estimator preserved theoretical accuracy for additive noise amplitudes up to 50% of the field amplitude, whereas the direct variance estimator showed substantial noise-induced bias and required correction. Bi-layer simulations and phantom experiments reproduced the predicted direction and onset of TOF-dependent decorrelation-rate trends in layered media. In vivo, the recovered blood-flow index tracked the expected TOF-dependent cuff-occlusion and reactive-hyperemia response in the forearm. During the single-subject Sudoku task, the left-forehead recording showed a TOF-dependent relative blood-flow-index increase of +0.8 {+/-} 1.9% at TOF = 400 ps, +9.8 {+/-} 2.2% at TOF = 600 ps, and +15.2 {+/-} 5.6% at TOF = 800 ps. This pattern is consistent with increased sensitivity to deeper tissue at longer photon path lengths, but requires cohort-level validation before quantitative interpretation as cognitive activation. Conclusions: Coupling temporal speckle-contrast analysis with swept-source iNIRS yields a proof-of-concept, depth-resolved platform for blood-flow sensing. By estimating TOF-resolved speckle contrast through the g_1-derived {kappa}^2 route, TOF-iSCOS suppresses additive-noise bias while preserving sensitivity to deeper dynamic tissue layers. The present single-channel results bridge continuous-wave iSCOS, interferometric NIRS and time-domain diffuse correlation spectroscopy (TD-DCS), and motivate future multi-channel and cohort studies for scalable cortical hemodynamic monitoring.

Optical fluorescence microscopy enables visualization of biological structures and dynamics. However, the intrinsic diffraction limit, especially axially, and depth-related scattering noise compromise the image resolution and fidelity. Computational 3D deconvolution is a promising approach for mitigating these issues, yet its execution is hindered by inaccurate and cumbersome theoretical modeling or experimental measurement of 3D point spread function (PSF), as well as ineffective 3D noise regularization. Furthermore, in the 3D super-resolution regime, there remains a lack of standardized tools for evaluating 3D super-resolution fidelity. Here, we present the 3D adaptive deconvolution and evaluation (3D-ADE) toolkit, which comprises 3D-Ada deconvolution with physics-oriented automatic 3D-PSF calibration, and 3D-SQUIRREL for 3D super-resolution quality assessment. It effectively resolves noise instability, eliminates the need for 3D-PSF calibration, and reliably assesses the fidelity of 3D resolution extension via deconvolution, physical, and deep-learning-based methods. Accessible via multiple software platforms, 3D-ADE enhances the versatility of 3D deconvolution and fills the gap in 3D super-resolution evaluation tools, and thereby advances volumetric fluorescence imaging applications.

Three-dimensional fluorescence microscopy often exhibits anisotropic resolution because axial information is poorly sampled and more blurred than lateral information, which complicates quantitative interpretation of fine 3D structures. Although optical remedies and computational restoration have been explored, many approaches require demanding system calibration or rely on accurate PSF models and assumptions that are difficult to satisfy across all samples and modalities. Here we present DeepIso, a self-supervised isotropy restoration framework that couples supervised pretraining with an internal-learning inference stage to estimate degradation directly from the measured volume. Without explicit PSF specification or enforced lateral-axial structural equivalence, DeepIso recovers axial frequency content and improves the continuity of elongated structures while retaining fine features, with superior performance over existing computational approaches in terms of both visual inspection and quantitative metrics. The method is validated on synthetic benchmarks and experimental datasets, demonstrating isotropy enhancement across confocal, light-sheet, and 3D structured illumination microscopy, thereby supporting downstream volumetric analysis including segmentation and tracking.

For decades, the photon counting histogram (PCH) was used as the sole method to quantify fluorophore numbers in a diffraction-limited focal volume. This technique combines spatial excitation profiles, and the distribution of photon counts to register the photon emission statistics of individual fluorophores. However, this approach has not yet been transferred to widefield fluorescent imaging due to the lack of fast and single photon sensitive camera sensors which can capture the photon emission statistics of a single fluorophore. Here, we explore avenues towards quantitative analysis of the active fluorophore number by leveraging recent advancements in single photon avalanche diode (SPAD) array technology. Binary exposures of a SPAD array can be synchronized with picosecond laser pulses to measure the PCH in a widefield setting. Then, by modeling the statistical relationship between the active fluorophore number and the PCH in a region of interest following a laser pulse, we can perform Bayesian inference of this number. The model is demonstrated experimentally by counting quantum dots and various numbers of fluorescent dye molecules bound to DNA origamis. We find that this method has several important applications in widefield microscopy, including enhanced localization microscopy and constrained fitting of multiple unresolvable fluorescent emitters.

We present a user-friendly correlative second harmonic generation (SHG) and atomic force microscopy (AFM) imaging workflow for quantifying the nanomechanical properties of collagen in unfixed, unlabeled tissue sections. SHG Aligned Profiling for Elasticity and Segmentation or SHAPES utilizes SHG imaging to guide AFM force mapping, enabling label-free, anatomically specific selection of regions of interest, facilitating spatially resolved characterization of fibrous collagen morphology and local stiffness. Compared to brightfield or confocal contrast-guided AFM mapping, SHG improves anatomical specificity without fixation or staining, enabling downstream analysis on the same tissue section while preserving spatial correspondence between structural (SHG) and mechanical (AFM) data. Rapid identification of regions of interest also reduces the risk of costly AFM tip breakage, improving throughput and reducing operator burden. Utilizing only standard turn-key commercial systems common in many user facilities and clinical laboratories, multimodal image coregistration and automated identification of anatomical regions of interest are employed to integrate SHG and AFM datasets across complex sample topographies. Coregistration of SHG and AFM images substantially increases the number of usable datasets per measurement session, facilitating translation of these complementary modalities and bridging nanomechanical imaging with clinical research practice.

This study presents the use of sub-micron nanobubbles (NBs) as contrast agents for ultrasound localization microscopy (ULM), a super-resolution imaging technique that visualizes microvascular structure and flow beyond the acoustic diffraction limit. While ULM has traditionally relied on micron-sized microbubbles (MBs), the reduced dimensions and prolonged circulation times of NBs make them attractive candidates for localization-based imaging. However, their weaker acoustic responses present significant challenges for reliable detection and tracking. To address this challenge, we developed the ULM Master GUI, an interactive framework for optimization of the complete ULM processing pipeline. Using custom ultrasound-compatible wall-less gelatin flow phantoms containing vessel-mimicking channels and bifurcations ranging from 100 to 500 m, we demonstrate that NB-based ULM achieves velocity reconstruction and flow partitioning measurements comparable to conventional MB-based ULM. Across all investigated geometries, NBs faithfully reproduced the underlying flow patterns and hemodynamic behavior despite their substantially reduced acoustic scattering. These findings establish the feasibility of NB-based ULM, expand the range of contrast agents available for localization microscopy, and provide a foundation for future super-resolution ultrasound imaging using nanoscale acoustic contrast agents. The ULM processing GUI is publicly available at https://github.com/grisha1998/ulm-super-resolution-toolbox.

Optical microscopy is fundamental to modern life-science research, yet interpreting its results requires precise modelling of point spread functions (PSFs) within complex environments. This manuscript introduces a versatile and efficient approach to wave-optical PSF calculations that extends existing frameworks by incorporating detection PSF modelling through the principle of reciprocity. Accompanying this work is a free MATLAB software package centred on a single, minimalistic core function, PlaneWaveExc.m, which utilizes a plane-wave superposition based on the Richards-Wolf model. Despite its simplicity, the framework accounts for "real-life" complexities such as systemic aberrations, arbitrary amplitude and phase modulations, and stratified media with complex-valued refractive indices. We demonstrate the softwares broad applicability through diverse case studies, including single-molecule imaging, STED microscopy, the segmented aperture of the James Webb Space Telescope, and coherent wide-field iSCAT microscopy. Each example is supported by dedicated scripts to facilitate adaptation for specific research needs.